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Abstract 


This  paper  reviews  the  newly  emerging  theory  of  a priori  error  bounds 
for  aggregation  of  mathematical  programming  models . The  main  technical 
building  blocks  are  described  in  some  detail,  and  the  available  applica- 
tions are  briefly  reviewed. 

V 


A.  AGGREGATION  THEORY  AND  ITS  APPLICATION  TO  MODELING 
IN  MATHEMATICAL  PROGRAMMING 

Aggregation  is  very  commonly  used  when  designing  mathematical  pro- 
gramming models.  One  is  ten^jted  to  say  that  nearly  all  models  are  aggregated 
in  one  way  or  another,  as  by  combining  individual  manufactured  items  into 
product  grotjqps,  employees  into  manpower  classes,  machines  into  machine  cen- 
ters, customers  into  customer  zones,  time  into  time  periods,  and  so  on. 

There  are  at  least  two  major  reasons  for  the  ubiquitous  use  of  aggrega- 
tion: it  simplifies  the  data  requirements  of  model -building,  and  it  simpli- 
fies the  task  of  solving  the  model  once  it  is  built.  Data  development  is 
simplified  because  fewer  data  elements  are  needed  and  because  "aggregate" 
data  are  more  amenable  to  development  using  statistical  sampling  methods  and 
summary  information  already  available  within. the  orgainization.  Aggregated 
models  are  easier  to  solve  because  they  are  smaller,  which  reduces  memory  re- 
quirements and  computing  time  , and  they  may  also  have  a more  tractable  mathe- 
matical structure  owing  to  the  disappearance  of  complicating  detail. 

It  is  surprising  that  so  important  and  comnonly  used  a modeling  device  as 
aggregation  is  almost  always  employed  in  an  ad  hoc  fashion  without  any  theor- 
etical justification.  Educated  intuition  and  common  sense  are  only  rarely 
assisted  by  analytical  argument.  It  would  almost  seem  that  the  old  saying 
"modeling  is  an  art  and  not  a science"  has  frightened  off  any  serious  study 
of  aggregation  theory  in  modeling.  Or  perhaps  the  technical  difficulty  of 
aggregation  theory  in  the  neighboring  fields  of  economics  and  statistics  [10] , 
and  the  apparent  inapplicability  of  that  theory  to  modeling  in  mathematical 
programming,  has  been  discouraging.  In  any  case,  ^ul  effort  is  long  overdue 
to  help  fill  the  void. 

We  wish  to  stress  that  our  concern  is  with  the  modeling  rather  thzm  algor- 


ithmic uses  of  aggregation  in  mathematical  programming.  The  latter  topic  has 
been  treated  in  depth  by  Lee  [11]  and  Zipkin  [12] . 
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what  is  required  is  a way  to  make  quantitative  error  estimates  for  al- 
ternative aggregated  model  designs  before  the  model  is  actually  built  and  run 
numerically.  These  error  estimates  would  supplement  common  sense  and  should 
enable  the  design  of  models  which  constitute  a provably  good  tradeoff  between 
accuracy  and  simplicity. 

More  specifically,  we  need  the  following  type  of  result.  Let  (P)  be  a 
fully  detailed  model  (before  cmy  aggregation) : 

(P)  Maximize  f(x,y)  subject  to  (x,y)  eF 

x,y 

where  x is  a vector  of  variables  not  to  be  aggregated  and  y is  a vector  of 
variables  to  be  aggregated  into  a new  vector  z.  Let  (Q)  be  an  aggregated 
version  of  (P) : 

(Q)  Maximize  f(x,z)  subject  to  (x,z)  eF 

x,z 

where  the  objective  function  f and  the  feasible  region  F have  been  altered 
to  accomodate  the  aggregation.  We  seek  priori  bounds  e^  and  such  that 

(1)  - 1.  v(P)  - v(Q)  , 

where  v(*)  denotes  the  optimal  value  of  an  optimization  problem.  By  "4 
priori",  we  mean  that  and  can  be  calculated  without  having  to  optimize 
(Q) . We  also  seek  a solution  recovery  mapping  Y(0  from  F into  the  space  of 
y such  that  if  (x,2)  is  5-  optimal-^  in  (Q)  and  yeY(x,z),  then  (x,y)  is 
(6+ej^+E2) -optimal  in  (P) . One  might  think  of  y(*)  as  a "disaggregation 
rule" . 

All  of  the  aggregation  theoretic  results  to  be  described  subsequently 
are  of  this  general  type. 

■^That  is,  C(x,2)  ^ v(Q)  - 6 . Clearly  0-optimality  corresponds  to  the  usual 
notion  of  (true)  optimality. 


4 


- 3 - 

A little  reflection  v5>on  (1)  and  the  role  of  the  solution  recovery 
mapping  should  convince  the  reader  that  they  are  exactly  what  is  required 
for  making  practical  decisions  about  which  aggregated  version  of  (P)  to 
adopt,  and  then  how  to  relate  aggregate  solutions  back  to  (P)  once  the 
aggregate  model  is  built  and  solved^  Notice  that  roost  per- 

tinent single  measure  of  aggregation  accuracy  — for  instance,  it  is  the 
width  of  the  interval  based  on  v(Q)  within  which  v(P)  must  lie.  Notice 
also  that  nothing  whatever  is  said  about  how  closely  a solution  to  (Q) 
approximates  a solution  to  (P)  in  "decision  space";  all  assertions  are  in 
"payoff  space".  We  believe  decision  space  proximity  to  be  much  less  relevant 
an  issue  in  a practical  sense  than  payoff  space  proximity. 

The  model  design  process  typically  involves  not  one  but  several  or  even 
many  possible  choices  for  (Q)  . Let  ^2^®^  denote  the  k priori 

error  bounds  as  a function  of  (Q) . Ideally,  we  would  like  the  following  to 
be  true  for  these  functions: 

a)  available  as  explicit  closed  form  expressions  without  the 
need  for  iterative  numerical  calculation, 

b)  equal  zero  in  the  special  cases  where  the  true  error  is 
known  to  vanish, 

c)  continuous  functions  of  the  data  of  (P) . 

Properties  b)  and  c)  together  imply  that  the  bounds  are  "close"  to  zero  if 
(P)  is  "close"  to  any  special  case  where  the  true  error  is  known  to  be  zero. 
We  would  also  like  the  solution  recovery  mapping  to  be  as  tractable  as 
possible . 

These  ideals  are  actually  achieved  for  most  of  the  applications  men- 
tioned later. 

^We  shall  not  address  here  the  issue  of  how  to  improve  the  ^ priori  bounds 
using  information  available  after  (Q)  has  been  optimized,  although 
usually  this  is  easy  to  do. 


Our  efforts  to  develop  this  type  of  aggregation  theory  for  mathematical 
progrcimming  have, up  to  the  present  time,  revolved  largely  about  specific 
applications.  These  applications  will  be  reviewed  (Sec.  6) , but  the  main 
body  of  this  paper  is  devoted  to  a review  of  the  four  analytical  techniques 
which  have  enabled  these  applications. 

1.  FIRST  ANALYTICAL  TOOL;  THE  RESTRICTIVE  APPROXIMATION  THEOREM 

Consider  a general  problem 
(P^)  Maximize  f (x)  subject  to  xeX 

and  the  following  restrictive  approximation  to  it 
(Q^)  Maximize  f(x)  subject  to  xeX  fl  X , 

where  X and  X are  both  subsets  of  the  same  decision  universe. 

Restrictive  Approximation  Theorem 

Assume  that  X is  not  en^Jty,  that 

(2)  f(x)  ^ f(x)  for  all  x in  X H x 

holds,  and  that  a Feasibility  Recovery  Rule  r(*)  exists  which  associates 
to  every  point  x in  X some  point  r(x)  in  X H x in  such  a manner  that 

(3)  f(r(x))  ^ f(x)  - e for  all  x in  X 
holds  for  some  0 . Then 

(4)  0 <_  v(P^)  - v(Q^)  ^ e 

and, for  any  6^0,  every  6-optimal  solution  of  (Q^)  will  necessarily  be  (6+e) 
-optimal  in  (P^) . 
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Proof . We  have 

v(Q^)  ^ Si;5>reinuin  f(r(x))  ^ v(P^)  - e , 
xeX 

where  the  first  inequality  follows  from  the  fact  that  r(x)  e X P x and  the 
second  inequality  follows  from  (3) . This,  with  the  evident  implication  of 
(2)  that  v(Q^)  ^ v(P^) , proves  (4) . It  remains  to  show  that  any  6-optimal 
solution  of  (Q^) , say  x°,  is  (6+e) -optimal  in  (P^) , that  is, 

(5)  f(x°)  > v(P^)  - (6+e). 

This  follows  from 

f(x°)  ^ f(x°)  ^ v(Q^)  - i5  ^ v(P^)  - E - 6 , 

where  the  first  inequality  follows  from  (2) , the  second  from  the  definition 
of  x°,  and  the  third  from  (4).  Q.E.D. 

This  result  slightly  generalizes  a theorem  first  proven  in  [7] . 

2.  SECOND  ANALYTICAL  TOOL;  lAGRANGEAN  DUALITY 
Linear  programming  duality  theory,  its  generalization  to  convex  pro- 
gramming [3]  and  its  formal  application  to  integer  programming  [4] , provide 
bounds  of  use  in  aggregation  theory.  The  possibilities  are  illustrated  with 
an  excunple  adapted  from  Zipkin  [12]. 

Consider  a general  linear  programming  problem 

2 t t t 

(P  ) Maximize  c x + d y subject  to  Ax  + By  ^ b 

x,y>o 

where  x and  y are  coliomn- vector  variables  and  the  data  arrays  c,d,A,B  and  b 
are  conformable.  The  y vector  is  to  be  aggregated  into  a scalar  variable  z 
according  to 

(6)  y = 6z>  where  0 > o is  a fixed  column-vector. 
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Upon  appending  (6)  to  (P  ) and  eliminating  y,  one  obtains  the  aggregated 


problem 


Maximize  c^x  + d^0z  subject  to  Ax  + B9z  < b^, 
x,z^O 


which  we  shall  presume  is  not  infeasible. 


Proportional  Restriction  Theorem  for  Linear  Programs 


Assume  that  some  bound  K is  known  such  that 


(7)  Zy*  <_  K 

j ^ 

2 

holds  for  at  least  one  optimal  solution  (x*,y*)  of  (P  ) . Corresponding  to 
each  variable  y^  , let  be  the  optimal  value  of  the  trivial  (continuous 
knapsack)  problem 


Maximize  d.-  u^B . subject  to  u^B6  > d^0 
u>o  3 => 


Define 


£ A Max  { 0 ,K  Mcix  { a . } } 

j ' 


O ^ v(P^)  - v(Q^)  ^ e 


Moreover,  if  (x°,z°)  is  6-optimal  in  (Q^)  for  some  6 ^0,  then  (x°,y®=  6z°) 

2 

will  necessarily  be  (6+e) -optimal  in  (P  ) . 


Proof. 


v(P^) 


2 

= v(P  l^jYj  hypothesis  (7) 


= Infimian 
u>o 


Supremum  c^x+d^y+u^  (b-Ax-By)  subject  to  EyAK 
x>o,y>o  j ^ 


by  duality  theory  [3] 


I 
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^ Supremum  c x+d  y+u  (b-Ax-^)  subject  to  Zy . ^ K 
x>o,y>o  j ^ 


where  u is  an  optimal  dual  solution  of  (Q  ) 


= u^b  + Supremum  (c^-u^A)x  + Supremum  (d^-u^)y  sub  j.  to  ^y 


x>o 


y>o 


; r 


Now 


u\>  = v(Q^)  by  the  definition  of  u and  LP  duality  theory,  the 


supremum 


over  X equals  O because  u must  satisfy  u^A  > c^  , eind  the  supremum  over  y is 


Max  {O ,K  Max  { d . - u B . } } 

j ' ' 


because  it  amounts  to  a simple  continuous  knapsack  problem. 
Thus  we  have  shown 


2 2 — t 

v(P  ) ^ v(Q  ) + Max  { O ,K  Max  {d . - u B . }}  . 

j ^ ^ 

But  a.  > d,  - u^.  by  (8)  and  the  fact  that 
j — j j 

“t  t “ 

u B6  > d 9 and  u > o 


hold  by  the  definition  of  u . Hence 


2 2 

v(P  ) < v(Q  ) + Max  {O  ,K  Max  {a.}}  , 

j ' 

which  proves  the  right-hand  inequality  of  (10) . The  left-hand  inequality  of 

2 2 
(10)  holds  because  (Q  ) is  really  just  a restriction  of  (P  ) . 

Suppose  now  that  (x°,z°)  is  6-optimal  in  (Q^)  for  some  6 ^ o . Then 
v(Q^)  - 6 <_  c^x°+  d^9z°  = c^x°+  d^y°. 


Adding  to  this  the  relation 


v(P^)  < e+  v(Q^) 


8 


from  (10) , one  obtains 


,„2,  . t o ,t  o 

v(P  )-6^E  + cx+dy  . 

This  result  shows  that  (x°,y°)  is  (6+e) -optimal  in  (P^) , since  (x°,y°) 


is  clearly  feasible  in  (P  ) 


Q.E.D. 


This  theorem  is  one  of  a family  of  similar  results  which  can  be  proved 
in  a similar  fashion.  For  example,  there  can  be  several  simultaneous  aggre- 
tations  of  the  type  (6)  on  different  subsets  of  variables,  there  can  be 
assumptions  different  from  (7)  concerning  y*,  and  (6)  could  be  replaced  by 
a more  general  linear  relation  in  which  6 is  a nonnegative  matrix  and  z is 
a vector. 

3.  THIRD  ANALYTICAL  TOOL:  THE  OBJECTIVE  FUNCTION  APPROXIMATION  THEOREM 


Consider  the  general  problems 


Maximize  f(x)  subject  to  xeX 


Maximize  f(x)  subject  to  xeX, 


where  X is  an  arbitrary  non-empty  set  and  f and  f are  both  real-valued 


functions  bounded  cibove  on 


The  function  f is  interpreted  as  an 


approximation  to  the  "real"  objective  function  f. 


Objective  Function  Approximation  Theorem  [ 5] 


Let  Ej^  and  be  scalars  (not  necessarily  nonnegative)  satisfying 


- E < f(x)  - f(x)  ^ E„  for  all  xeX. 


It  is  permissible  for  f and  f to  assume  the  value  - so  long  as  they 
both  do  it  on  the  same  proper  subset  of  X. 


Then 


I 


I 


» 


> 


(12)  ■ - ^2 

and,  for  any  6^0/  6-optimal  solution  of  (t^)  . i ■• 

(6+ej^+e2)  “Optimal  in  (P^)  . 

4.  FOURTH  ANALYTICAL  TOOL;  PROJLCTI  Jt. 

Consider  the  general  problem 

4 

(P  ) Meucimize  f(w,y)  subject  to  (w,y)eG 

w,y  weW,yey. 

The  "projection"  of  this  problem  "onto  the  space  of  ti.-  • ■ • 

defined  as : 

4 ^ 

(Q  ) Maximize  ^(w)  subject  to  weW ' F 

w 

where 

(13)  ^(w)  A Supremum  f(w,y)  subject  to  (w,y) 

Y 

(14)  F A any  svibset  of  { w:  the  problem  defi.ni.ii^  < - 

feasible  solution-. 

We  adopt  the  convention  that  ^ (w)  = - “>  for  those  values  of  • ' . • 

defining  problem  in  (13)  is  infeasible.  Note  that  the  oondi'-.  t ■ 

4 

implicit  feasibility  constraint  which  must  hold  for  (P  ) as  . • 
Projection  Theorem  [2] 

If  (w°,y°)  is  E-optimal  in  (P^) , where  e ^0  , then  w°  auat  i- 
in  (q"^)  . If  w°  is  e '-optimal  in  (Q^)  and  y°  is  e"-opti»al  in  rr..  - 

defining  (zi(w°)  , where  e'  ^0  and  e"  ^0  > then  (w°,y°)  is  (i  '*  “ - 

4 4 4 4 

(P  ) . If  problem  (P  ) is  feasible,  both  (P  ) and  (0  ) have  ident.ca  • • 


rental  values. 


10 


♦ 


> 


I 


I 


» 


I 


I 


It  can  also  be  shown  [2]  that  projection  is  a convexity-preserving 

4 

transformation.  That  is,  if  (P  ) has  a concave  objective  function  over  a 

4 

convex  feasible  region,  then  the  same  will  be  true  of  (Q  ) . 

Projection  is  closely  related  to  the  notion  of  hierarchical  decisions 
and  control  (e.g.,  see  [9]),  although  previously  published  work  along 
these  lines  has  not  fully  exploited  this  fact. 

One  of  the  most  valuable  uses  of  projection  is  as  a device  to  facili- 
tate the  application  of  each  of  the  three  previously  described  devices. 

This  is  illustrated  in  the  next  section. 

5.  OBJECTIVE  FUNCTION  APPROXIMATION  CUM  PROJECTION 

The  four  analytical  tools  are  building  blocks  which  can  be  combined  in 
various  ways  to  obtain  various  aggregation  theoretic  results.  In  this  sec- 
tion we  show  how  projection  cind  objective  function  approximation  cam  be 
combined  in  a fruitful  manner. 

Start  with  a feasible  problem  of  the  form 

(P^)  Maximize  f(x)+g(y)  subject  to  (x,y)eW,  xeX,  yeY 

x,y 

into  which  some  new  aggregation  variables  z are  introduced  along  with  an 
associated  set  of  definitional  constraints  z = Z(x,y) : 

(P  ) Maximize  f(x)+g(y)  subject  to  (x,y)eW,  xeX,  yeY,  z = z(x,y) . 

x,y,z 

For  instamce,  each  component  of  z might  be  defined  as  the  sum  of  selected 
con?)onent8  of  y.  Since  the  aim  is  to  replace  y by  the  aggregate  variables 
z,  project  (P^)  onto  the  space  of  the  x and  z variables: 

7 

(P  ) Maximize  f(x)+  ^(x,z)  subject  to  xcX  and  (x,z)eF 

x,z 


I 


where 


(15)  iHx,z)  A Supremuro  g(y)  subject  to  (x,y)eW,  yeY,  z=»  Z(x,y) 

y 

(16)  F A a suitable  subset  of  {(x,z)  : the  problem  defining  0(x,z) 

has  a feasible  solution}. 

By  a "suitable”  subset  in  (16)  we  mean  F should  be  chosen,  if  possible,  so 

as  not  to  spoil  whatever  special  structure  (P^)  may  have  that  could  be  of 

7 

use  in  solving  the  aggregate  problem.  Finally,  let  us  replace  in  (P  ) by 
some  more  tractable  approximation  to  obtain 

(Q^)  Maximize  f(x)+  fS(x,z)  subject  to  xeX  and  (x,z)  eF. 

x,z 

What  is  the  relationship  between  the  original  problem  (P^)  and  the 
aggregated  problem  (Q^)?  The  answer  can  be  found  by  examining  the  successive 
steps  of  the  transformation 

(P^)  ^ (P®)  (p'^)  (Q^) 

in  light  of  the  results  established  in  the  preceding  sections.  The  relation- 

5 6 5 

ship  between  (P  ) and  (P  ) is  obvious,  although  technically  (P  ) could  be 

6 6 7 

viewed  as  a projection  of  (P  ) . The  relation  between  (P  ) and  (P  ) is 

available  from  the  Projection  Theorem  (the  role  of  w is  played  by  (x,z)). 

7 5 

The  relation  between  (P  ) andCQ  ) is  available  from  the  Objective  Function 
Approximation  Theorem.  Putting  it  all  together,  we  obtain  the  following 
end  result. 

Theorem 

Suppose  that  f+ff  is  bounded  above  on  the  set  of  feasible  solutions  of 
(Q^)  , and  that  f+fS  equals  - ® if  and  only  if  f+0  equals  - * on  a (possibly 
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empty)  proper  subset  of  this  set.  Suppose  further  that  scalars  and 
exist  such  that 

- t^<_  i»;x,z)  - (z((x,z)  ^ 

holds  for  all  (x,z)  such  that  0 is  finite  and  feasible  in  (Q^) . Then 

(17)  - <_  v(P^)  - v(Q^)  ^ • 

Moreover,  if  (x°,z°)  is  5-optimal  in  (Q^)  for  some  6^0  and  y°GY  satisfies 
o o 

(x  ,y  ) ew  and 

(18)  g(y°)  ^ 0(x°,z°)  - , 

then  (x°,y°)  is  (6+g^+  G2)-optimal  in  (P^) . 

7 

Proof . The  Objective  Function  Approximation  Theorem  applied  to  (P  ) and 
(Q^)  yields 

(19)  ■ ^1  - - ^2  ’ 

5 7 

But  v(P  ) = v(P  ) by  the  Projection  Theorem,  so  (17)  is  at  hand.  Suppose 
that  (x°,z°)  is  6-optimal  in  (Q^)  and  that  y°GY  satisfies  (x°,y°)  gW  and 
(18) . Then  (x°,y°)  is  feasible  in  (P^)  and 

f(x°)+g(y°)  ^ f (x°)+iJ<(x°,z°)  - g^^  by  (18) 

^ v(Q^)  - 6-  by  the  6-optimality  of  (x°,z°)  in  (Q^) 
^ v(P^)  - <5  “ by  (17) 

= v(P^)  - (6+ej^+  e^) » 

which  shows  that  (x°,y°)  is  (6+  Gj^+  G2)“ophimal  in  (P^)  . Q.E.D. 
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6.  APPLICATIONS 

Most  applications  to  date  have  been  in  the  area  of  distribution  system 
planning  emd  related  models,  but  recently  results  have  also  been  achieved 
for  a production  scheduling  model  and  a production/marketing  model. 

In  [7]  we  applied  the  Restrictive  Approximation  Theorem  to  achieve  ll 
priori  error  bounds  for  customer  aggregation  in  a class  of  distribution 
system  design  models.  Customer  aggregation  is  universally  practiced  in  this 
context  because  most  firms  have  many  tens  of  thousands  of  customers;  to  in- 
corporate all  of  them  individually  would  lead  to  excessive  data  development 
costs,  inordinate  computing  time  and  primary  storage  requirements,  and  an 
excessive  degree  of  solution  detail.  Our  error  bounds  deal  with  the  most 
natural  type  of  customer  aggregation,  namely  pro  rata  by  demand,  for  single 
commodity  one  or  two-stage  distribution  systems  with  possible  economies- 
of-scale  in  facility  throughput  costs.  A compeuiion  paper  [8]  extends 
these  results  to  the  multicommodity  case  with  the  help  of  an  auxiliary 
model  based  on  Lagrangean  relaxation.  It  turns  out  that  commodity 
aggregation  error  bounds  are  also  obtained  as  a byproduct  of  the  analysis. 

Zipkin  [12]  was  able  to  rederive  the  ik  priori  bound  of  [7] , for  the 
special  case  of  the  classical  transportation  problem,  using  the  approach  of 
Sec.  2.  He  also  presents  a priori  and  2k  posteriori  error  bounds  for  more 
general  network  models  aggregated  in  an  imaginative  variety  of  ways. 

Another  paper  [6]  applies  Projection  and  the  Objective  Function 
Approximation  Theorem  to  obtain  A priori  error  bounds  for  a class  of  pro- 
curement problems  in  a fairly  general  logistical  setting.  The  items  being 
aggregated  were  the  various  commodities  being  procured.  An  interesting 
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aspect  of  the  results  obtained  is  that  linear  progranming  can  be  used  to  de- 
termine the  best  aggregated  model  once  the  items  to  be  aggregated  are 
specified,  where  "best"  means  that  values  of  certain  arbitrary  coefficients 
in  the  model  are  determined  so  as  to  yield  the  smallest  a priori  error 
bound.  This  is  apparently  the  first  instance  where  LP  has  been  used  to  de- 
termine an  optimal  aggregative  approximation  to  a class  of  mathematical 
programming  problems . 

In  some  unpublished  notes  we  have  applied  Projection  eind  the  Restrictive 
Approximation  Theorem  to  obtain  a priori  error  bounds  for  product  aggrega- 
tion in  a machine  loading  and  scheduling  model  now  being  ir5)lemented  for  ein 
injection  molding  coit5>any.  The  model  involves  several  thousand  products 
made  using  several  hundred  tools  (mold  bases)  mounted  on  a hundred  machines. 

The  aim  is  to  decide  how  much  of  each  product  to  make  and  to  sell  during 
each  time  period  over  a full  season,  cind  also  to  decide  how  best  to  utilize 
tools  and  machines.  There  are  constraints  on  tool  eind  machine  availability, 
the  maximum  production  rate,  and  minimum  and  maximum  sales.  The  costs  for 
machine  operation,  inventory,  and  sales  are  linear  but  production  costs  are 
concave . A priori  error  bounds  and  solution  recovery  results  enable  the 
problem  size  to  be  reduced  while  preserving  the  special  structure  on  which 
the  specialized  solution  technique  depends  (this  technique  uses  Lagrcuigean 
releucation  [4]  in  conjunction  with  an  efficient  network  flow  optimization 
code  [1] ) . 

In  ^ulother  set  of  unpublished  notes  we  applied  the  Restrictive 
Approximation  Theorem  to  regional  aggregation  in  a comprehensive  planning  model 
now  being  implemented  by  a large  agricultural  supply  cooperative.  The  model 

I 

must  determine  how  much  of  each  conmodity  to  produce  where,  including  trades  and 
exchanges  with  con^titors,  and  how  great  a market  share  to  obtain  in  each 
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marketing  region  via  controllable  market  development  expenditures.  There 
are  pervasive  nonlinearities  in  the  market  development  cost  functions,  but 
otherwise  the  model  is  of  linear  programming  type  (more  precisely,  it  could 
be  viewed  as  a generalized  network  with  side  constraints) . Our  aggregation 
theoretic  results  apply  to  the  selective  coarsening  of  the  original  grid  of 
customer  regions,  which  was  taken  for  convenience  as  a checkerboard  of 
square  plots  50  miles  on  a side. 

7,  (CONCLUSION 

The  results  achieved  to  date  show  clearly  that  a priori  error  bounds 
for  aggregation  in  mathematical  programming  models  is  a topic  which  is  both 
tractable  mathematically  and  applicable  to  problems  of  practical  interest. 
Perhaps  the  greatest  need  at  this  time  is  for  empirical  studies  aimed  at 
assessing  the  quality  of  the  a priori  bounds.  Favorable  indications  here 
would  justify  the  development  of  additional  theoretical  tools  and  an  effort 
to  popularize  aggregation  theory  as  a model  design  technique. 

On  the  basis  of  the  preliminary  computational  evidence  currently  avail- 
able, it  would  appear  that  a priori  bounds  often  show  that  negligible  error 
is  associated  with  quite  substantial  aggregations.  Frequently  these  aggre- 
gations are  ones  an  experienced  modeler  would  have  accepted  on  intuitive 
grounds  anyway,  but  he  should  nevertheless  be  pleased  to  have  a solid  ration 
ale  for  use  when  defending  the  model's  assumptions.  As  the  extent  of  aggre- 
gation becomes  still  greater,  the  empirical  error  may  for  a time  remain 
acceptably  small  even  though  the  h.  priori  bound  becomes  unacceptable.  Effec 
tive  model  design  in  this  range  awaits  the  development  of  priori  error 
estimators  aimed  at  probedsle  range  rather  than  worst  case. 
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